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In a previous publication, a two-dimensional Green's function has been derived for a 
circular cylinder whose surface impedance around the periphery deviates from a constant 
value by a sinusoidal variation of small amplitude a. Here, this solution is evaluated 
asymptotically in the illuminated region under the assumption that the cylinder radius is 
large compan^d with the wavelength of the incident field. The asymptotic result is inter- 
preted in terms of geometrical optics generalized to apply to cylindrically curved, convex 
reflection gratings, and comprises the first and higher-order diffracted rays associated with 
a plane grating, together with geometrical divergence coefficients accounting for the surface 
curvature. General properties of the spectrum of reflected rays are observed. The behavior 
of the first-order diffracted rays, and in particular their domain of existence as a function of 
surface periodicity, is discussed in detail. 



1. Introduction 

In the theory of diifraetiou at high frequencies, a 
lowest-order approximation to the scattered field 
solution can be constructed by the methods of 
geometrical optics. If a plane wave represented by 
a continuum of parallel rays falls on a perfectly con- 
ducting, gently curved convex object, the reflected 
geometrical optics field comprises rays which emerge 
from the obstacle surface at the angle of incidence; 
a reflected ray lies in the plane of incidence and has 
an amplitude along the ray direction governed by 
a divergence coefficient which accounts for the 
spreading, due to the surface curvature, of a narrow 
reflected ray bundle. If the surface properties are 
characterized by a constant surface impedance Z^, 
the geometrical optics approximation assumes that 
the reflected ray amplitude is, in addition, deter- 
mined by a reflection coefficient which is the same as 
for an infinite plane surface with impedance Z^ tan- 
gent to the object at the refiection point. To 
verify the validity of these approximations for a sur- 
face with constant curvature, one may examine the 
rigorous solution to the problem of diffraction by a 
constant impedance cylinder in the high-frequency 
limit [Wait, 1959]. If the surface properties of the 
scatterer are not constant but are characterized by a 
spatially varying surface impedance Z^, geometrical 
optics predicts that the scattered field can be con- 
structed as before provided that one now employs 
the appropriate impedance at the specular reflection 
point in calculating tlie ray reflection coefficient 
[Keller, 1956]. Implicit in all geometric-optical 
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approximations is the requirement that any varia- 
tions in surface characteristics — whether changes in 
the radii of curvature, or in the local impedance 
properties — arc ^ 'small'' within an interval of length 
equal to the wavelength of the incident field. 

To gain a deeper insight into the simultaneous 
influence of surface curvature and variability of 
surface properties on the field diffracted by an object 
at high frequencies, the authors have studied the 
problem of diffraction by a large circular cylinder 
whose surface impedance around the periphery 
deviates from a constant value by a sinusoidal varia- 
tion of small amplitude a (see (2)). Apart from the 
relative simplicity of the mathematical analysis lead- 
ing to a solution of this problem, the sinusoidal 
variation is highly desirable since an appropriate 
choice of the periodicity permits the simulation of 
either slow monotonic changes, or of rapid fluctua- 
tions. In the former range, a conventional geo- 
metric-optical interpretation of the dominant portion 
of the high-frequency diffracted field is expected to 
apply. Its domain of validity as a function of both 
the surface curvature and the rapidity of impedance 
variation can be assessed by studying the high-fre- 
quency asymptotic solution for arbitrary values of 
periodicity. In this manner, one may arrive at 
estimates for the magnitudes of the ''small" terms 
mentioned at the end of the preceding paragraph, 
which delimit the range of applicability of geo- 
metrical optics. 

The above discussion suggests a division of the 
asymptotic analysis of the rigorous solution into two 
parts : (a) the evaluation of the high-frequency fields 
for arbitrary periodicity of the surface impedance, 
and (b) the transformation of this solution into the 
previously mentioned geometric-optical form when 
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the periodicity is very slow. Part (a) of this 
program is carried out in the present paper wherein 
it is shown that the asymptotic solution of the 
scattered field is interpretable via the geometrical 
optics of a curved, convex reflection grating, and 
that it can be constructed in terms of the specular 
and diffracted rays familiar from the theory of the 
plane grating, modified by geometrical divergence 
coefficients to account for the surface curvature 
(sec. 2). An interesting feature is the dependence 
of the domain of existence of each diffracted ray on 
both the periodicity of the surface impedance and on 
the radius of curvature, as discussed in section 4. 
Phase (b) of the investigation, presented in the 
companion paper, is concerned with the case where 
the period of the surface impedance variation is very 
large compared to the incident wavelength, thereby 
obscuring the periodic aspect of the surface proper- 
ties. As noted above, one then expects the reflected 
field to be derivable from geometric-optical consider- 
ations applied to surfaces with locally constant 
characteristics. The reducibility of the first-order 
asymptotic solution to the geometric-optical formula 
deduced from a ^ 'local reflection'^ argument can thus 
be employed to furnish some bounds on the range 
of applicability of conventional geometrical optics 
when used in connection with gently curved surfaces 
possessing a variable surface impedance. 

The formal solution to the diffraction problem 
mentioned above [Felsen and Marcinkowski, 1962, 
and Marcinkowski and Felsen, 1961. These refer- 
ences will henceforth be called I and II, respectively. 
A related study for a slightly noncircular cylinder 
has been carried out by Clemmow and Weston, 1961] 
has been obtained by a Green's function procedure 
wherein the problem of determining the amplitudes 
in a harmonic series expansion of the scattered field 
Gs is reduced to that of solving an inhomogeneous, 
second-order difference equation with variable co- 
efficients. The solution of the latter has been found 
by assuming the expansibility of Gs as a power series 
in the small perturbation parameter a. While we 
discuss in detail only the first-order perturbation 
(sec. 3 and appendix), the form of the asymptotic 
solution is readily ascertained to any order in a and 
leads to the expressions in section 2. Stress is 
placed througliout on the physical interpretation, 
rather than the detailed derivation, of the mathe- 
matical results; the interested reader will find 
additional information in reference II. 

The physical configuration and choice of coordi- 
nates are shown in figure 1. While the circular 
cylindrical coordinates (p, </>, z) are employed in the 
derivation of the formal solution and its asymptotic 
representation, the coordinates Bf, ^„ and s are 
convenient for the geometric-optical interpretation 
of the result. A plane wave of wavelength X with a 
magnetic field H^ along the positive z axis is incident 
on the cylinder along the ray OQ identified by the 
angle ^\ The field reflected by the cylinder is 
observed in the (p, <^) plane at the point P{p, 0). 
On the surface of the cylinder (p=a), the magnetic 
field Hz(p, 0, 00 is required to satisfy the boundary 
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Figure 1. The geonielric optic parameters. 
condition 



1 dH,{p,M') 

ioie dp 



= -Z(4>), 



(1) 



with the variable surface impedance Z(<^) given by 



Z(<t>) =Z [^ J' [1 +a COS />(<^-^o) ] • 



(2) 



e and jLt are the permittivity and permeability of the 
medium, and Z is an arbitrary complex number. 
For^a passive surface impedance, the restriction 
Re Z(0)>O is imposed. In (2), a small, sinusoidal 
variation with complex amplitude a (0<|q:K<1) is 
superimposed upon the constant impedance 
Z(/z/€)^/2 Ti^g periodicity of the surface properties 
is described by a finite, positive integer ]) = 2TralL, 
where L is the spatial period of the impedance 
variation. This provides the useful parameter 
p/ka=\/L involving the propagation constant 
k=2'7r/\. For the sake of convenience the angle of 
observation and the arbitrary real phase angle 
00 are restricted to the ranges — 7r<C(0— 0O<C7r and 
— TT < (00—00 < TT. The units are rationalized MKS, 
and a time dependence exp(— ^co^), where i is the 
imaginary unit, co the angular frequency, and t the 
time, is suppressed. As mentioned previously, the 
formal solution of this boundary value problem has 
been presented in references I and II. We proceed 
now to its asymptotic evaluation in the illuminated 
region of the cylinder. 

2. Spectrum of Reflected Grating Rays 

The total axial magnetic field H^ in the illuminated 
region may be decomposed into incident and reflected 
contributions. 



G{p, 0, cf>')=Gi(p, Cj), Ct>') + Gr{p, Ct>, 00 , 



(3) 
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where the subscripts i and r identify the incident and 
reflected parts, respectively. The incident plane 
wave is given by 

Gi{p, <i>, ct>') = e-'''^ ^"^ ^^-'^'\ 0<6/,<7r/2. (4) 

A detailed asymptotic evaluation of Gr by the 
saddle point method has been carried out to 0(a) 
and is summarized in the appendix and in section 
3. However, the main properties of the asymptotic 
solution to any order in a can be inferred from the 
first-order results by repeated application of a re- 
cursion relation derived in references I and II. 
These more general relations are presented in this 
section. Their interpretation in geometric-optical 
terms is facilitated by the use of the geometric- 
optical coordinates 6iy dry s in figure 1. 

6i is the angle between the incident ray M^M and 
the normal ON to the cylindrical surface at the 
point of reflection M(0<^i<7r/2). dr is the angle 
between the reflected ray MP and the same normal 
ON (— 7r/2<^,<7r/2). The coordmate s is the 
straight-line distance along a reflected ray MP from 
the point of reflection M to the point of observation P. 
To permit a simple interpretation of the asymptotic 
results, it is necessary to admit nonspecularly re- 
flected rays so that the angle of reflection dr need 
not be equal to the angle of incidence di] this gen- 
eralization of conventional geometrical optics allows 
the inclusion of the diffracted rays associated with a 
periodic structure or grating, dr is considered posi- 
tive when the reflected ray lies on the side of the 
normal ON opposite the incident ray as shown in 
figure 1; when both incident and reflected rays lie 
on the same side of OA^ then Br is considered negative. 
If 6r<i—Bi and s is sufficiently large, the observation 
point P{py (j)) may lie below the line OQ. Under 
these conditions the ray MP is reflected from a 
point M on the cylinder where 0— 0'>O and is 
observed at a point P in space where </>— </)'<^0. 

With these definitions, the reflected field in the 
illuminated region may be shown to comprise a 
spectrum of reflected waves determined as follows: 



Grip, <t>y <!>')= S Ao(np)^r.{p)D{np)e''^^'^^^ (5) 
where 

^n{p) = J^ R(n, m,p) a^^+'^i, m-0, 1,2,... 



D(np)=[l+'M1 
r(np)—a 



-1/2 



(np)J 
cos drjnp) 

ddiinp) 



(6) 



and 



sin 6r{np) = sm 6i{np) ±n jy 6r,^0) (7) 



with < ^, inp) <7r/2 and -7r/2 <^, (np) <7r/2 . 



Each reflected ray associated with each integer 
n will in general emerge from a different point M 
in figure 1 [with a different di{np)] to reach the same, 
prescribed observation point P. The exponent in 
the term Ao(np) gives the phase of the incident 
plane wave at each reflection point, while i^n{p) 
gives the reflection coefficient associated with each 
reflected ray. It is well known that a parallel 
bundle of rays incident on a cylinder produces a 
divergent bundle of reflected rays [Bremmer, 1949 
and Keller, 1956]. These rays appear to come from 
a point located by the focal length r(np) which is 
measured along that part of the segment MP in 
figure 1 extended behind the reflection point 
M {r(np)=MF). The resulting decrease in ampli- 
tude along a reflected ray is described by the diver- 
gence coefficient D(np). D(np) and r(np) above 
are generalizations of analogous geometric optical 
parameters for the rays reflected specularly from a 
cylindrical surface with constant impedance. In 
tlie derivation of these parameters from purely 
geometrical considerations, the cylindrical surface 
and its tangent have been assumed to be con- 
tinuous functions of p and 0, and Br{np) and its 
first derivative have both been assumed to be 
continuous functions of diinp). The term ks(np) 
in the exponential gives the phase of a reflected 
ray propagating from the point of reflection M to 
the observation point P(p, </>). The saddle point 
evaluations show that the angles of incidence and 
reflection for an nth order reflected ray are related 
by the well-known grating law in (7), familiar from 
the theory of a plane grating. The angle 6^ in (7) 
is the angle NOQ which the normal ON in figure 1 
makes with the incident ray (l)^(\dr,\ = di(np), ^,>0 
in fig. 1). The explicit form of Aoinp), D{np) and 
r{np) in (5) and (6) is obtained easily by comparing 
the integral expressions for the higher-order (in a) 
reflected contributions (see ref. II) with the 0{a) 
expressions in the appendix. For the explicit 
evaluation of the reflection coefficients ^nip), how- 
ever, a detailed asymptotic evaluation of the higher 
order integral contributions is required and is not 
included here. 

From (5) it is evident that a finite number of 
reflected rays reaches a prescribed point P(p, </>). 
The orders of the reflected rays range from — % to 
+^2, where % and n2 are nonnegative integers which, 
in general, are not equal. They are the largest 
values of n for which (7) can be satisfied subject 
to the indicated restrictions on di(np) and 6r{np). 
The perturbation procedure produces reflection 
coefficients ^nip) which are expressed as power 
series in a involving terms of the form a^'^+l'^'. The 
coefficients of these series, and therefore the reflec- 
tion coefficients .^nip), are independent of the 
radius of the cylinder. Therefore, it is reasonable 
to expect that these reflection coefficients will be 
identical with those obtained for a plane surface 
with the same sinusoidally modulated surface im- 
pedance, where the plane surface is taken to be 
tangential to the cylinder at the point of reflection. 
This supposition has been verified to O(a^) by com- 
paring our reflection coefficient in (12) with that 
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obtained by Hessel [1960] for the plane, sinusoidally 
modulated impedance surface. For an nth order 
reflected ray (n>0), the leading term of the reflec- 
tion coefficient ^n(p) is of Oia"^) in the perturbation 
parameter. For any prescribed 0(q:'^) there are a 
total of 71+1 rays of 0(a'^), a result which may 
be derived from the recursion relation appearing in 
the solution (ref. I). 

It is significant to emphasize that the asymptotic 
reflected field representation in (5), expressed in 
terms of the geometric-optical coordinates (Ot, dr, s), 
depends only on the local properties of the surface 
in the vicinity of the point of reflection of a given 
ray. Hence it appears reasonable to assume that 
the dominant contribution to the reflected field in 
the illuminated region of a general, gently curved 
cylindrical convex surface with periodic impedance 
properties can be constructed as in (5), provided 
that the constant cylinder radius ^'a^\ appearing in the 
focal length r{np) and hence in the divergence coeffi- 
cient D(np), is replaced by the radius of curvature of 
the surface at the point of reflection of a given ray. 
For curved surfaces with nonperiodic, slowly varying 
properties, the above assumption was put forth by 
Keller [1956] and verified by Keller and others 
[Keller, Lewis, and Seckler, 1956] for various spe- 
cific obstacle configurations, in which the reflected 
field comprises the usual specularly reflected rays 
only. The results herein suggest a direct extension 
to include as well the nonspecularly reflected rays 
associated with periodic surface properties. 

3. Reflected Rays to Order a^ 

To O(a^) in the perturbation, the reflected fields 
in the illuminated region consist of three rays whose 
properties are described here: 

+a[Gl{p,^,ct>' ,p) + Gl{p,4>,cl>' -p)]. (8) 

In terms of the geometric optic parameters, the 
reflected field to 0{a°) is given asymptotically by: 



with a reflection coefficient 

cos di—Z 



Uo= 



cos di-\-Z 



o<e,<^, (9) 



(10) 



For this specularly reflected ray (n=0), the simpler 
notation di{np)-^di is used, together with analogous 
simplifications for the other relevant parameters: 
dr = di and ddrfddi = 1. Any point P(p, </>) in the 
lit region may be reached by a specularly reflected 
ray for which 0<^^<7^/2. These results are, of 
course, well-known from the theory of scattering 
by a constant impedance cylinder [Imai, 1954; 
Franz, 1957, 1954; Franz and Beckmann, 1956; 
Keller, 1956]. 

The corresponding solution of O(a^) in (8) com- 
prises two reflected waves for which n = ±1. The 



field for 71 = +1 is given asymptotically by (see 
appendix) 

Gl(p,ct>,<t>',p)=A,(p)E,(p)D{p)e''^^^\ (11) 

where the reflection coefficient is found to be 



R.(P) 



_ —Z cos gi(y) g'"'*"" ~'^''' 

~[Z+cos eAp)][z+cos eTipV] 



(12) 



The angle at the point of reflection is <^ (p) = 
(t>^±di(p) for <9^$0. Each of the two reflected rays 
of O(a^) in (8) satisfies the grating law given by (7), 
with n=±l. From the grating law it follows^ that 
the ray geometry of these two reflected rays is 
invariant to the simultaneous substitutions p^>—p 
and dr,-^—dr,. As a result, the ray geometry is 
symmetrical about the incident ray </>' in figure 1 
regardless of the phase angle </)o in (2) provided that 
a ray with the index n= + l(—l) for ^,;>0 is com- 
pared with a ray of index n= — l( + l) for 6rj<C0, for 
the same angle of incidence di(p) = \drj\. However, 
the symmetry of the reflected fields is, in general, 
different from the symmetry of the corresponding 
ray geometry due to the presence of the reflection 
coefficient Ri(p), given to O(a^) by (12). This 
coefficient depends upon the phase angle 0o through 

the term exp[i^(0(^)— </)o)] = exp[ip (</>'— 0o±^z(^))] 
for ^^^0. If 00 has the special values ^(^' — 
0o)=±^7r, where m=0 or a positive finite integer, 
then both the ray geometry and the reflected fields 
are symmetrical about the incident ray <^'. This 
symmetry property follows directly from the 
symmetry of the variable surface impedance pre- 
scribed by (2). To O(a^) in the perturbation and 
asymptotically in ka ior ka >>1, it is easy to verify 
that (8) satisfies the variable impedance boundary 
condition in (1). 

If the periodicity parameter p is small or if the 
angle of incidence 6i(p) is large, the angle of reflec- 
tion for both grating rays in (7) will be positive 
(dr(±p) >0) . With respect to the surface normal ON 
in figure 1, reflected rays of this kind lie on the side 
opposite the incident ray. This situation will occur 
simultaneously for both rays n=±l in (7). This is 
illustrated in Regions A and D of figure 2, where the 
three reflected rays to O(a^) in the perturbation have 
been drawn in four different situations which can 
arise in the lit region. On the other hand, if the 
periodicity parameter p is large or if the angle of 
incidence di(p) is small, the situation may arise 
where Ori—pXO for ^,>0. This is iflustrated in 
figure 2 for ray 3 in Region B and ray 2 in Region C. 
Let us consider the behavior of ray 3 as the surface 
normal ON passing through the point of reflection 
moves from Region A to B. The angle of reflection 
^r{—p) changes continuously from a positive value in 
Region A to a negative value in Region B. When 
the surface normal ON lies on the dashed line which 
bounds Regions A and B, then dr{—p)='~0. Under 
this special condition the reflected ray lies along the 
surface normal ON. Similar conclusions hold for 
ray 2 in Regions C and D. From (7), the angle of 
incidence 6i(p)=do appropriate to the special 
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ON -SURFACE NORMAL 



sine„= -^ .9,-20' 



Fir, I IRE 2. Rejlected rays to order a. 
limiting condition dr{ ±p) =0 is given by the equation 



L 



(13) 



Figure 2 lias been drawn to scale for ^o = 20°. 

Let us consider an angle of incidence such that 
all three reflected rays are in Region A. As the 
normal ON in figure 2 is varied continuously from 
Region A to D, all three rays change their position 
continuously. It is evident that the ray geometry 
is symmetric about the incident ray 0' providing the 
roles of rays 2 and 3 in figure 2 are reversed upon 
reflecting about <^^ If :p=0, the three rays coalesce 
into one specularly reflected ray (dr=di) and Regions 
B and C disappear. For this special condition, a 
single reflected ray is obtained from (8) by adding 
the three reflected rays into one specularly reflected 
ray. To 0{a^) it is easy to verify that the resulting 
reflected ray gives the expected asymptotic solution 
for te>>l for a constant impedance cylinder with 
the surface impedance Z/ = Z(fx/ey^^{l-^a). 

4. Domains of Existence of the Reflected 
Rays of Order a^ 

The domains of existence of the reflected grating 
rays are those illuminated or ^'lit'^ regions wherein 
we find geometric-optical field contributions com- 
prising real rays. For these regions, real angles 
of incidence and reflection exist which simultaneously 
satisfy the grating law given by (7), and the associ- 
ated limitations on 6i{np) and dr(np). In view of 
these limitations, the lit regions for the grating 
rays will depend upon \/L and may be considerably 
smaller than the lit region for specularly reflected 
rays. For the two reflected rays of order a^ this 





(a) -^=0 



i^)T = T 





(c) ^ = 1 



REFLECTED RAY n = -l 



REFLECTED RAY n = + l 



(d)r = T 



-LIMITING RAY n = + l 
■LIMITING RAY n = -l 



P=^ >0 



= 2 



Figure 3. The domains of existence of the real 
reflected grating rays n=ibl. 

dependence is examined in this section as a function 
of \/L, the relative periodicity of the variable surface 
impedance. The corresponding properties of the 
higher order grating rays m.ay be obtained by 
appropriate generalizations of these simpler results. 

Eacli of tlie two grating rays n= ± 1 has its own 
lit region. As described previously in connection 
with figure 2, the region illuminated by the ray 
n= + l is the mirror image of the region illuminated 
by the ray 7^= — 1 reflected about the incident ray 
</)'. The ranges of the lit regions associated with 
the two reflected rays n=±l are shown in figure 3 
for \/L=0, 1/2, 3/2, and 2. (It should be noted 
that the saddle point evaluations, and therefore the 
asymptotic solutions herein, are not valid on the 
limiting rays or in the transition regions surrounding 
the limiting rays.) 

For the constant impedance cylinder, \/L=0. 
For this special condition, both rays n=±l have 
the same lit region as the specularly reflected ray 
71=0 (fig. 3a). As \/L increases from \/L=0 to 
\/L=2, the extent of the lit regions associated with 
each ray progressively decreases as shown in figures 
3a to 3e. When \/L=2, each of the two lit regions 
degenerates into a single line as shown in figure 3e 
and the associated reflection coefficient for each ray 
vanishes from (12). If X/i>2, the grating law 
cannot be satisfied for real angles, and the lit regions 
for the rays n= ±1 disappear. 

For each ray family 7^= ±1 there are two charac- 
teristically different limiting rays which define the 
shadow-lit boundaries. For one limiting ray, the 
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incident ray always arrives at grazing incidence, 
di(:hp) = Tr/2. Since the angle of incidence is fixed, 
the corresponding reflected ray always emerges from 
the same point on the cylinder with an angle of 
reflection which varies with X/i in the manner 
prescribed by the grating law. There is one limiting 
reflected ray of this type for each of the two ray 
families n=±l as shown in figure 3. For the other 
limiting ray, the roles of the angles of incidence and 
reflection are reversed and the reflected ray always 
leaves the cylinder surface at a grazing angle (grazing 
reflection) . In this instance the angle of incidence, 
and therefore the point of reflection, varies with 
X/L according to the grating law as shown in figure 
3. There is one limiting ray of this type for each 
ray family n=ztl. 

The domains of existence of the reflected rays 
n=±l described above can be understood from a 
consideration of the first-order diffracted rays on 
a plane reflection grating which are governed by the 
grating law in (7). In fact, the various regions 
depicted in figure 3 are precisely those which would 
be obtained for a plane grating of the same relative 
periodicity \/L if the incident wave sweeps out the 
entire range from normal to grazing incidence. 
Consider, for example, the case \/L=S/2 in figure 
3d whence for the ray n= — l, with ^„>0, sin dr= 
sin di—3/2. The diffracted ray does not emerge from 
the plane of the grating until sin^^>l/2; for the 
maximum angle of incidence (i.e., grazing incidence) 
sin di=l, whence sin ^,= — 1/2. Thus, the domain 
of existence of the first order diffracted ray is the 
angular region 7r/2>— ^;.>7r/6 (see fig. 4), with 
—dr=7r/2 and —dr=T/6 corresponding to di = 7r/6 
and di=Tj2, respectively. This illuminated region 
for the plane grating is exactly equivalent to that 
shown for the ray n= — l in figure 3d, if the curved 
periodic surface is approximated at the point of 
reflection by a plane surface. This is in accord with 
the previously noted quasi-optic expectation (ka 
»1) that the field at any point in the illuminated 
region of a curved surface with periodic properties 
can be constructed from the geometric-optic results 
for a plane surface wherein diffracted rays arising 
due to the periodicity are included. The effect of 
surface curvature is accounted for by a divergence 
coefficient. This discussion of the domains of 
existence of the illuminated regions has been restrict- 
ed to the two grating rays arising in the perturbation 
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FiGUEE 4. The plane grating for X/L=3/2. 



solution to O(a^). Nevertheless, the conclusions 
obtained here may be easily extended to the higher 
order rays by making the substitution \/L^n\/L 
for an nth order ray. Consequently, if \/L=l/2, 
then figure 3 may also be interpreted as representing 
the extent of the illuminated regions for the specular 
rays and for the first four orders of the grating rays 
(ji=l, 2, 3, 4). In this case, no real diffracted rays 
exist for n^4:. The regions outside those shown 
illuminated in figure 3 represent '^grating shadow'' 
regions whose properties remain to be investigated. 

5. Appendix 

This appendix contains a summary ^ of the deri- 
vation of the geometric-optical formula (11) from the 
rigorous integral expression for the reflected wave 
Gl, given by (see ref. II) 



Z ^p 

Gl{p,4>,<l>',p)=^e 



where 



X 






H/,kp) 



a{v)'a{v+p) 
a(v) = H',{ka) +iZH,(ka) . 



dv (14) 



Hv indicates a Hankel function of the first kind and 
the prime denotes differentiation with respect to the 
argument. If the appropriate large argument 
asymptotic approximations are substituted for the 
Hankel functions, then this integral becomes 



where 



G\{p,ct>,<i>' ,V)^ ^ Q{v,v) e'^^^'-^' dv (15) 



f{v,p)={^'—<j>-~~'^) ^^^ 7+sin 7 



— (sin 7i— Ti cos 7i) (sin 72—72 cos 72) (16) 

P P 



Q{v,v)= 



ip( 0'-0o — M-h^- 



(Z+sin7i)(Z+sin72) 



f sin 7i sin 72 "]^/^ . . 
L27r^psin7j ^ ^ 



V v+v 



C0S7i= 



ka 



0<Re (7, 7i or 72)<7r. 



(18) 



The contour c runs from — 00 to + 00 passing 
above the real v axis for ^Q{v)<^—ka and below for 
Re(j/)>/:a. The integrand in (15) contains a rap- 
idly varying exponential term and the slowly vary- 
ing function Q(v,p). Therefore, it is in a suitable 
form for an asymptotic evaluation by the saddle 
point method when /:p>/:a>> 1 (for simplicity it 

2 Details of the evaluation can be found in reference II. 
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is assumed that arg(/:) = 0). 
given by (18) it follows that 



Fr 



cos 72 = COS 7i + 



ka 



the definitions 



(19) 



The contour c hi (15) may be deformed into a 
steepest descent path passing through the saddle 
point v = vq on the real v axis. The equation for the 
ocation of this saddle point is 



-0— ^— 7o+7io+T2o="0 



(20) 



where 70, 710, and 720 are defined in (18) by letting 
The asymptotic approximation to G\ is thus given 

by 



0\' 






Q{v,,p)e'^^^^'^'^' (21) 



which hivolves the angles 70, 710 and 720. Ap- 
propriate transformations (whose details are omitted 
here) replace (18) by the simpler restrictions 0<^i<^ 
x/2, — 7r/2^^,<^7r/2 defining the illunhnated region, 
and furthermore 



ikpj{vQ,p)^iks{p) —ika cos Bi-\-ip I - ± di \ 

where 

s{p)-==-p sin Jo— a sin 710 

= pcos[±{cl>'-ct>)+d,+dr] 



(22) 



(23) 



When evaluated at the saddle point, the relation 
(19) transforms into the grathig law given by (7). 
The focal length r previously defined is given here 
by the equation 



r{p)=a- 



sm 7io 



cos drip) 



1- 



sm 720 
sin 7io 



1 + 



(16, 



(24) 



while the divergence coefficient D arises from 



^(z»=,-^ 



■\ 1/2 



-sin 7io 



sin Yo sin 720 V (>'\f {v,P) ~\ 



(25) 



where 



miv.v) '] ^1/1 



7o ^«in7io tt sin 720, 



.(26) 



By means of (16) to (26) the as3nnptotic I'esult in 
(21) may be transformed into the simple geometric 
optical formula (11). The approximate expressions 
used in (15) for the Hankel functions break dow^n in 
the transition regions where di or |^^|«7r/2. To de- 
scribe the transition phenomena in the narrow re- 
gions surrounding the limiting ra3"s in figure 3, it is 
necessary to employ a more detailed analysis than 
that presented above. 
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